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Abstract. Fermat Last Theorem, which inspired mathematicians during 300 years, 
is proved by Andrew Wiles. Even among mathematicians there is a narrow circle 
of specialists, who can read this proof and understand all details. Is it a reason 
for pessimism? No, since arithmetics if entire numbers contains broad variety of 
problems with a simple statement, which might be not less intricate. One of them 
arises in elementary geometry. 



1. Elementary problem on pyramid. 

Primary education (Elementary School) and secondary education (High School) 
in Russia are united into one stage that now lasts 11 years (from 6 year old to 17 
years old). Mathematics is among disciplines studied during these years. Below we 
consider a problem, which can be suggested to 10-th or 11-th year students in the 
course of geometry. It is typical, though is a little more complicated than usual. 

Problem on pyramid. In triangular 
pyramid ABCD three sides of triangle 
ABC in its base are given: 

\BC\=a, \CA\=h, \AB\^c. 

From corners A and B two perpendiculars 
are drown to the faces BCD and ACD 
respectively. Their lengths are given: 

\AF\ = /, \BC\ = g. 

Find the length of the segment [FG] con- 
necting feet of these two perpendiculars 
[AF] and [BC] . 

Let's consider in brief the steps leading 
to the solution of this problem. First we 
draw all three heights in triangle ABC: 
these are segments [^-ff], [Sif], [CA/]. It's known that they cross at one point. 
Denote it by L. For the sake of simplicity we consider the case when triangle 
ABC is acute-angled. In this case point L lies inside the triangle ABC . Now let's 
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apply Pythagor's theorem to rectangular triangles AMC and BMC. As a result 
we obtain the system of equations with respect to the length of the segment [AM] : 

( \AM\^ + \CM\^ = \AC\'^, 
^^■^^ { {\AB\-\AM\f + \CM\'^ = \BC\'^. 

Solving the system of equations (1.1), for lengths of [AM] and [BM] we get 

\ABl+jACl-jBCl 
l^^l - 2\AB\ ' 

^^■^^ \AB\^ + \BC\^ - \AC\^ 

l^^l = 2\AB\ • 

Similar formulas can be obtained for \AK\, \KC\, \BH\, and \HC\: 



\AK\ 
\CK\ 
\BH\ 
\CH\ 



\AC\^ + \AB\^ - 




2\AC\ 




\AC\^ + \BC\^- 




2\AC\ 




\BC\- + \AB\- - 


\ACf 


2\BC\ 




\BC\^ + \ACf - 




2\BC\ 



Let's replace \AB\ — \AM\ by \BM\ in the second equation of the system (1.1). 
Then we can derive the following formula for the length of segment [CM]: 



(1.3) \CM\ = 



\AC\^ + \BC\^ - \AM\^ - \BM 
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Similar formulas can be derived for the lengths of segments [AH] and [BK] : 



\AH\ = 

(1.4) 

\BK\ = 



|AB|2 + |BC|2 - \AK\^ - \CK\'^ 



In order to calculate lengths of segments [KL] and [HL] we use similarity of trian- 
gles: AKLC ~ AM AC and AHLC ~ AMBC. This yields: 



(1.5) 
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Now let's draw segments [FH] and [Gif]. According to the theorem on three 
pcrpcndicnlars, we have GK _L AC and FH _L BC. Then, since we already know 
\AF\ and \BG\, we can calculate lengths of segments [FH] and [GK]: 

[FH[ = y^[AC\^ - \CH\^ - \AF\^, 

[GK\ = ^/\BCf^^]CK^^^]BG^. 

In order to derive first of these two expressions we applied Pythagor's theorem to 
rectangular triangles AHC and AHF. Second expression is derived by Pjdihagor's 
theorem applied to triangles BKC and BKG. 

Orthogonal projections of the points F and G onto the plane of the base of 
pyramid belong to the straight lines AH and BK. Denote these projections by F 
and G respectively. For the sake of simplicity we consider the case when points 
F and G are above the base of pyramid (i. e. in upper halfspace separated by the 
plane ABC), and when their projections F and G belong to the segments [HL] 
and [KL] respectively (see Fig. 1.2 and Fig. 1.3). Due to similarity of triangles 



l\GKG ~ /\KBG and /\HFF ~ /\FAH we derive the following formulas: 
(1.6) ]GG[ = ^t^, WF[ = 1^"^" "^^1 



[BK\ ' ' ' [AH] ' 

(1.7) \KG] - ]HF\ - 

The length of the segment [Gi^] (see Fig. 1.4 below) is determined by cosine theorem 
applied to the triangle GLF: 



(1.8) ]GF\ = |iG|2 + |iF|2 - 2 ]LG[ \LF\ cos{KLH) 

Note that angles ZKLH and ZKCH complete each other to a straight angle. 

Indeed, triangle LKC is rectangular (sec Fig. 1.4 below). Same is true for triangle 
LHC. Hence for the angles of these two triangles we can write the equalities: 

'KLC + laJL = 9Q°, 
HLC + HCL = 90°. 
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Adding these equalities and taking into account that KLC + HLC = KLH and 
'KCL + HCL = KCH, we get the required equahty KLH + KCH = 180°. As an 



immediate consequence of this equality we can write the equality for cosines: 

cos(KLH) = -cos{KCH). 

Cosine of the angle KCH can be determined by applying cosine theorem to the 
triangle ABC, which lies in the base of pyramid ABCD: 

= 2\AC\\BC\ • 

Lengths of segments [LF\ and [LG] in formula (1.8) can be calculated as follows: 

\LF\ = \LH\ - \HFl \LG\ = \LK\ - \KG\. 

This is obvious from Fig. 1.2 and Fig. 1.3. Now the length of segment {FG\, which 
was to be found, is calculated by Pythagor's theorem (see Fig. 1.5): 



\FG\ = sJ\FG\' + {\GG\-\FF\Y. 



So problem on pyramid is solved. This is typical stereometric problem that 

can be used to test the knowledge of some; basic fa(;ts and spatial imagination of 
students. Its solution considered just above is not tricky. But it is rather huge, and 
we cannot write simple explicit formula expressing \FG\ through parameters a, h, 
c. /, and g. Therefore we should give numeric values for these parameters, choosing 
them so that they provide simple numeric values for ultimate result and for results 
of all intermediate calculations. Thus, another problem arises, problem of choosing 
proper numeric values for a, b, c, /, and g. We shall consider this problem below. 
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2. Orthogonal matrices. 

Let's apply the coordinate method to the problem on pyramid. Here we have two 
natural triples of orthogonal vectors. First consists of vectors AF, FH, HC, second 
is formed by vectors BG, GK, and KG. Let's consider three unitary vectors ei, 
62, 63 directed along vectors AF, FH, and HG. Then choose other three unitary 
vectors directed along vectors BG, GK, and KG. Vectors ei, 82, 63 and hi, h2, h3 
form two bases consisting of unitary vectors orthogonal to each other. Such bases 
are called orthonormal bas6s (ONB). Let's consider the following expansions 
binding vectors of two ONB's: 



(2.1) 



hi = ^S'f6fc, 



SI 


Si 






Tl 


rpl 


rpl 

^3 


Sf 


Si 


Si 


T = 


n 




rp2 


Sf 


SI 


Si 






n 


n 



Coefficients of the expansions (2.1) are usually arranged into square matrices, which 
are called treinsition matric6s: 



(2.2) 



Matrices S and T in (2.2) implement direct and inverse transitions from base to 
base, they are inverse to each other, i. e. their product is a unitary matrix: 

S -T = T ■ S = E. 

If we treat (2.1) as transition from the base ei, 62, 63 to the base hi, h2, h3, then 5* 
is called dir6ct transition matrix, while T is called inv6rse transition matrix. 

Matrices S and T in our case are binding two ONB's. Therefore components of 
these matrices arc bound by a series of relationships. If we denote by 5* and T* 
transposed matrices, i. e. if we denote 





s\ 


si 


sf 




Tl 


rp2 


rp3 


S' = 


Si 


Si 


SI 


rpt ^ 


n 


rp2 

^2 


Ti 




SI 


Si 


Si 




n 


rp2 

^3 


Ti 



then these relationships for components of S and T can be written as follows: 
(2.3) S* ■S = E, T* ■T = E. 

From (2.3) and from S-T = T-S = Ewe immediately derive S* = T and T* = S. 

Matrices that satisfy the relationships (2.3) are called orthogonal matric6s. 
Sum of squares of elements in each column and in each string of orthogonal matrix 
is equal to 1. So we have the relationships 



(2.4) 



J2{Sl f = = 1 all ^ = 1' 2. 3. 
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Sums of products of elements from different columns and/or different string are 
equal to zero. This property is expressed by the relationships 

3 3 

(2.5) J2SISI = J2S^S1 = for M <?• 

i=l i=l 

The relationships (2.4) and (2.5) are easily derived from (2.3). Moreover, from (2.3) 
one can derive the following relationships for determinants of 5 and T: 

{detSf = l, (detr)2 = l. 

Therefore detS* = detT = ±1. Looking attentively at Fig. 1, one can note that 
AF, FH, HC and BG, GK, KG arc oppositely oriented triples of vectors: first is 
left, while second is right. Hence bases ei, 62, 63 and hi, h2, h3 are also oppositely 
oriented. This fact is reflected by the sign of determinants of transition matrices: 

(2.6) detS' = detr= -1. 

Further we shall be interested in the case when all components of matrices S 
and T are rational numbers. Components of S belonging to the same column can 
be brought to common denominator, and hence, they can be written as 

^'-dl' ^l-rfT' ^l-rfT' 

^^■^> ^^'T,' 

Ql _ ''I C2 _ '^2 r,3 _ rs 

From (2.4) for entire numbers pi, p2, Ps, and di in (2.7) we derive the relationship 

(2.8) {pif + {P2f + {psf = {d,f. 

If four entire numbers satisfy the relationship (2.8), we say that they form Pytha- 
gorean tetrad. Each column in orthogonal matrix with rational components is 
related with some Pythagorean tetrad of entire numbers. Thus, in (2.7) we have 
three Pythagorean tetrads determined by transition matrix S: 

(2.9) {pi,P2,P3,di), (91,92,93,^2), (r-i,r2,r3,rf3). 

Pythagorean tetrads of entire numbers (2.9) are orthogonal to each other in the 
sense of the following relationships: 

Pi 9i + P2 92 + Pa 93 = 0, 
Pi n + P2 r-2 + Pa rs = 0, 
n qi + r2 92 + r-3 93 = 0. 
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In order to determine an orthogonal matrix it's sufficient to have two orthog- 
onal Pythagorean tetrads, for instance, {pi,P2,P3,di) and (91, 92, 93, £^2)- Third 
Pythagorean tetrad then will be determined by the relationships 



d3 



P2 P3 
92 93 



d3 



Pi P3 

91 93_ 

di d2 



d3 



Pi P2 

91 92^ 

di d2 



This is the consequence of the fact that third vector in orthonormal bases (ONB's) 
are determined by vector product of first two vectors: 



63 = -[ei, 62], 



hs = [hi, ha]. 



The difference in sign here is due to the condition (2.6), which expresses difference 
in orientations of bases ei, 62, 63 and hi, h2, h^. 

Returning to the problem on pyramid, we use the fact that vectors AF, FH and 
HC are coUinear to base vectors ei, 62, 63: 

(2.10) AF = ai-ei, FH = a2 ■ e2, HC = as ■ 63- 
Similarly, vectors BG, GK, and KG are collinear to base vectors hi, h2, hs: 

(2.11) BG = Pi-hi, G^ = /32-h2, Kd = p3-h3. 
Vector BG is collinear to vector 63, while vector AG is collinear to vector h3: 



(2.12) 



BG = w • 63, 



AG = a- hs. 



Let's choose parameters uj and a in (2.12) to be rational numbers, and then let's 
apply the relationships (2.1). As a result for vector AG we obtain two expansions: 



(2.13) 



AG = AF + FH + HG, 

AC = a- {Si ei + S| 62 + 5| 63). 



Substituting (2.10) into (2.13) and comparing two expansions (2.13), we obtain 

(2.14) ai=aSl a2 = a Sl as = a Sl 

In a similar way, from (2.10) and (2.13) due to (2.1) and due to the expansion 
BG = BG + GK + KG we can derive the following three relationships: 



(2.15) 



(}2= OJT, 



3 ' 



P3=UJTI 



If components of transition matrix S are rational numbers, then components of 
inverse transition matrix T = are also rational. Therefore from (2.14) and 
(2.15) we obtain rationality of numeric coefficients a2, a2, as and /3i, 02, (is in 
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(2.10) and (2.11). This, in turn, provides rationality of lengths of segments [AC], 
[BC], [AK], [CK], [BH], [CH], [BG], [GK], [AF], and [FH]. 

Let's consider the vector FG, length of which is a final result in the problem on 
pyramid. For this vector we have an expansion: 

(2.16) FG = FH + HC -KG -GK. 

Let's substitute (2.10) and (2.11) into the expansion (2.16). This yields 

FG = a2 62 + as 63 - /?2 h2 - Ps hs = 62 + 03 63 - 

- 02 (^2 ei + Si 62 + Si 63) - 133 {Si ei + 5| 62 + 5| 63). 

Now it's easy to see that vector FG has rational coordinates in orthonormal base 
(ONB) formed by vectors ei, 62, 63. Therefore its length, in the worst case, is 
simplest irrational number obtained as a square root of rational number. The 
same is true for lengths of segments [AB], [AM], [BM], [CM], [KL], [HL], [AH], 
[BK], as well as for the lengths of segments [FF], [HF], [GG], and [KG]. For \AB\ 
this follows from the equality AB = Ad -Bd. Further we use formulas (1.2), 
(1.3), (1.5), (1.4); then formulas (1.6) and (1.7). Main conclusion that wc draw 
from what was said above is the following: orthogonal matrices with rational 
components give the algorithm for choosing numeric values of parameters a, b, c, 
f, g in the problem on pyramid such that wc get simple final result in this problem 
and simple results in all intermediate calculations. 

3. Constructing orthogonal matrices 

WITH rational components. 

Constructing orthogonal matrices with rational components is a separate prob- 
lem. First we consider regular algorithm for constructing such matrices. It is based 
on elementary rotations. Let's consider three entire numbers pi, P2, d, and suppose 
that they are bound by the relationship 

(3.1) (Pl)2 + (p2)2 = d\ 

Such numbers form Pythagorean triad. In contrast to Pythagorean tetrads they 

are well-known. There is a regular algorithm for constructing all Pythagorean 
triads (see, for instance, [1]). If r is a greatest common divisor of pi, P2, and q, 
then pi = T ■ pi, P2 = T ■ p2, d = T ■ d. From (3.1) we derive 

{Pif + {hf = d\ 

If Pi is even and p2 is odd, then d is odd. According to the regular algorithm 
described in [1], in this case we have the following expressions: 

Pi = 2 (rr? -\- m — r? — n), 
P2 = 4:mn + 2m + 2n + l, 
d = 2 {m^ + m + + n) + 1. 
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Here m and n are two arbitrary entire numbers. So Pythagorean triads arc parame- 
terized by three arbitrary entire numbers: m, n, and r. Suppose that we have some 
nonzero Pythagorean triad (pi,P2, d). Then we can consider two rational numbers 
Pi/d and P2/d, sum of their squares being equal to unity: 

Hence in half-open interval [0, 27r) there exists some angle (p such that 
(3.2) cos<^ = — , 

Angle (f in (3.2) is uniquely determined by numbers pi, P2, and d forming Pythago- 
rean triad {pi,P2,q)- Let's use this angle in order to define four matrices: 





cosip 


sin if 
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(3.3) 
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Matrices Sl^\ S\^\ S}^^ are geometrically interpreted as matrices of elementary 
rotations to the angle (p around coordinate axes. They arise as transition matrices 
in (2.1) in that case when base hi, h2, ha is got from base ei, 62, 63 by one of such 
elementary rotations. Matrix S* in (3.3) is interpreted as a matrix of inversion. It 
arises as transition matrix in (2.1) when 

hi = -ei, h2 = -62, ha = -63. 

All four matrices (3.3) are orthogonal. This can be checked by substituting 
them into (2.3). If (p is determined by the relationships (3.2), then all components 
of matrices (3.3) arc rational. Product of two orthogonal matrices is an orthogonal 
matrix (it's well-known that such matrices form a group). Therefore, choosing n 
Pythagorean triads and determining angles (^1, . . . , <p„ by relationships (3.2), we 
can consider the following product of corresponding matrices (3.3): 

n 

(3.4) s = {s*r ■ l[{si;^r ■ {s'^lf- ■ {s'^T- 

i=l 

Here e, a^, 7j are entire numbers either equal to zero or to unity. Formula (3.4) 
gives an algorithm for constructing orthogonal matrices with rational components 
in dimension 3. We shall call it a regular algorithm. 



P2 

svmp = —. 
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4. Some generalizations and open questions. 

Leonard Eulcr considered a class of matrices, which is a little more wide than 
class determined by the relationships (2.3). Euler's class is formed by matrices S 
with entire components that satisfy the following condition: 

(4.1) S'-S = N-E. 

Here N is some positive entire number. Matrices of Euler's class are called entire 
orthogonal matrices, number N is called a norm of orthogonality. If A?^ is a 

square of entire number, i. e. if A'' = M^, then matrix • 5 is an orthogonal 

matrix with rational components in the sense of standard definition by formulas 
(2.3). Leonard Euler has suggested an algorithm for constructing entire orthogonal 
matrices in the dimensions 3 and 4. His algorithm is described in book [2] . In papers 
[3 5] Euler's algorithm was generalized for n x n matrices in arbitrary dimension n. 
Due to the existence of two algorithms we have a scries of quite natural questions. 

— How do Euler's algorithm relate with regular algorithm, which is expressed by 
above formula (3.4) ? 

— Can we construct an arbitrary orthogonal matrix with rational components by 
Euler's algorithm ? 

— Is there the expansion (3.4) for an arbitrary orthogonal matrix with rational 

components, i. e. can it be constructed by regular algorithm? 

Answers to these questions are unknown to the author of this paper. Author will 
be grateful for any information concerning subject of this paper. 
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